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Unveiling systematic biases in ID LTE 

excitation-ionisation balance of Fe for FGK stars. A novel 
approach to determination of stellar parameters. 
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ABSTRACT 

We present a comprehensive analysis of different techniques available for the spec- 
troscopic analysis of FGK stars, and provide a recommended methodology which ef- 
ficiently estimates accurate stellar atmospheric parameters for large samples of stars. 
Our analysis includes a simultaneous equivalent width analysis of Fe I and Fe II spectral 
lines, and for the first time, utilises on-the-fly NLTE corrections of individual Fe I lines. 
We further investigate several temperature scales, finding that estimates from Balmer 
line measurements provide the most accurate effective temperatures at all metallicites. 
We apply our analysis to a large sample of both dwarf and giant stars selected from 
the RAVE survey. We then show that the difference between parameters determined 
by our method and that by standard ID LTE excitation-ionisation balance of Fe re- 
veals substantial systematic biases: up to 400 K in effective temperature, 1.0 dex in 
surface gravity, and 0.4 dex in metallicity for stars with [Fe/H] ~ —2.5. This has large 
implications for the study of the stellar populations in the Milky Way. 

Key words: stars: abundances — stars: late-type — stars: Population II 



1 INTRODUCTION 

The fundamental atmospheric (effective temperature, sur- 
face gravity, and metallicity) and physical (mass and age) 
parameters of stars provide the major observational foun- 
dation for chemo-dynamical studies of the Milky Way and 
other galaxies in the Local Group. With the dawn of 
large spectro scopic surveys to s tudy indiv idual stars, such 
as S EGUE jYannv et all I2009T). RAVE (|Steinmetz et all 
l200fj), Gaia-ESO (|Gilmore et all 120121 1. and HERMES 



Holmbcrg, Nordstrom, fe Andersen! 120071: lAn et all 120091: 
lArnadottir. Feltzing. fe Lundstrom 2010; Casagrande et al.l 



iBard en et al.|[2008l ). these parameters are used to infer the 
characteristics of different populations of stars that comprise 
the Milky Way. 

Stellar parameters determined by spectroscopic meth- 
ods are of a key importance. The only way to ac- 
curately measure metallicity is through spectroscopy, 
which thus underlies photometric calibrations (e.g., 



* Current address: Lund Observatory, Department of Astron- 
omy and Theoretical Physics, Box 43, SE-22100, Lund, Sweden: 
greg@astro.lu.se 



l201ll ), while high-resolution spectr oscopy is also used t o cor- 
rect the low-resolution results fe.g.. lCarollo et aTll2010h . The 
atmospheric parameters can all be estimated from a spec- 
trum in a consistent and efficient way. This also avoids the 
problem of reddening inherent in photometry since spectro- 
scopic parameters are not sensitive to reddening. The spec- 
troscopic parameters can then be used alone or in combi- 
nation with photometric information to fit individual stars 
to theoretical isochrones or evolutionary tracks to determine 
the stellar mass, age, and distance of a star. 

A common method for deriving the spectroscopic at- 
mospheric parameters is to use the information from Fel 
and Fen absorption lines under the assumption of hydro- 
static equilibrium (HE) and local thermodynamic equilib- 
rium (LTE). Many previous studies have used some vari- 
ation of this technique (e.g., ionisation or excitation equi- 
librium) to determine the stellar atmospheric parameters 
and abundances, and henceforth distances and kinematics, 
of FGK stars in the Milky Way. For example, some have used 



© 2012 RAS 



2 Ruchti et al. 



this procedure to estimate the effective te mperature, sur- 
face gravity, and metalticity of a sta r (e.g., iFulbright 2000; 
iProchaska et al.ll2000l ; iJohnsonll^Ool ) , whiie others use pho- 
tometric estimates of effective temperature in combination 
with the ionisation equiiibrium of the abundance of iron 
in LTE to estimate surface gravity and metallicity (e.g. , 
Mc William et al.ll 19951; iFrancois et alj|2003t iBai et alj|2004 



had signal-to-noise ratios greater than 100 : 1 per pixel. The 
equivalent widths (E Ws) of bo t h Fel and Fell lines , take n 
from the line lists of iFulbrightl l|2000r i a nd I Johnson! i|2002|) . 



were measured using the ARES code |Sousa et al.l 



2007) 



Allende Prieto et al.ll2006l ; lLai et alj|2008h . 

However, both observational (e.g.. iFuhrmannl 1998; 



llvans et al.ll200ll ; lRuchti et al.ll20l"ll;lBruntt et al.ll2012T) and 
theoretical evidence (e.g . , Thevenin fe Idiartl 19991 ; lAsplundl 
120051 ; iMashonkina et all 1201 if ) suggest that systematic bi- 
ases are present within such analyses due to the breakdown 
of the a ssu mption of LTE. More recently. [ Bergemann et al.l 
i|2012bf ) and lLind. Bergemann. fe Asplundl (|2012h quantified 
the effects of non-local thermodynamic equilibrium (NLTE) 
on the determination of surface gravity and metallicity, re- 
vealing very substantial systematic biases in the estimates 
at low metallicity and/or surface gravity. It is therefore ex- 
tremely important to develop sophisticated methods, which 
reconcile these effects in order to derive accurate spectro- 
scopic parameters. 

This is the first in a series of papers, in which we develop 
new, robust methods to determine the fundamental param- 
eters of FGK stars and then apply these techniques to large 
stellar samples to study the chemical and dynamical prop- 
erties of the different stellar populations of the Milky Way. 
In this work, we utilise the sample of stars selected from 
the RAVE survey originally published in lRuchti et al.l |201ll . 
hereafter Rll) to formulate the methodology to derive very 
accurate atmospheric parameters. We consider several tem- 
perature scales and show that the Balmer line method is the 
most reliable among the different methods presently avail- 
able. Further, we have developed the necessary tools to apply 
on-the-fly N LTE correction^ to Fe I lines, utilisi ng the grid 
described in lLind. Bergemann. fe Asplundl ^20 1 21). We verify 
our method using a sample of standard stars with interfero- 
metric estimates of effective temperature and/or Hipparcos 
parallaxes. We then perform a comprehensive comparison 
to standard ID, LTE techniques for the spectral analysis of 
stars, finding significant systematic biases. 



2 SAMPLE SELECTION AND OBSERVATIONS 

NLT E effects in iron are most prominent in low-metallicity 
stars (Lind. Bergemann. fe Asplundll2012l ; lBergemann et al.l 



12012b ). We therefore chose the metal-poor sample from Rll 
for our study. These stars were originally selected for high- 
resolution observations based on data obtained by the RAVE 
survey in order to study the metal-poor thick disk of the 
Milky Way. Spectral data for these stars were obtained us- 
ing high-resolution echelle spectrographs at several facilities 
around the world. 

Full details of the observations and data reduction of 
the spectra can be found in Rll. Briefly, all spectrographs 
delivered a resolving power greater than 30,000 and covered 
the full optical wavelength range. Further, nearly all spectra 



"NLTE correction" refers to the difference between the abun- 
dance of iron computed in LTE and NLTE obtained from a line 
with a given equivalent width. 



However, during measurement quality checks, we found that 
the continuum was poorly estimated for some lines. We 
therefore determined EWs for these affected lines using hand 
measurements. 



3 STELLAR PARAMETER ANALYSES 

We computed the stellar parameters for each star using two 
different methods. 

In the first method, which is commonly used in the liter- 
ature, we derived an effective temperature, Tlte-Fc, surface 
gravity, log (/lte-Pc, metallicity, [Fe/H]LTE-Fe, and microtur- 
bulence, wlte-Fo, from the ionisation and excitation equilib- 
rium of Fe in LTE. This is hereafter denoted as the LTE-Fe 
method. We used an i terative proc edure that utilised the 
M00G analysis program (ISnedenlll973T ) and ID, plane-parallel 
ATLAS-ODF model atmospheres from Kurucqj computed un- 
der the assumption of LTE and HE. In our procedure, the 
stellar effective temperature was set by minimising the mag- 
nitude of the slope of the relationship between the abun- 
dance of iron from Fel lines and the excitation potential of 
each line. Similarly, the microturbulent velocity was found 
by minimising the slope between the abundance of iron from 
Fe 1 lines and the reduced EW of each line. The surface grav- 
ity was then estimated by minimising the difference between 
the abundance of iron measured from Fel and Fell lines. 
Iterations continued until all of the criteria above were sat- 
isfied. Finally, [Fe/H]LTE-p c was chosen to equal the abun- 
dance of iron from the analysis. Our results for this method 
are described in Section [4] 

The second method, denoted as the NLTE-Opt method, 
consists of two parts. First, we determined the optimal effec- 
tive temperature estimate, Trjpt, for each star (see Section [5] 
for more details). Then, we utilised M00G to compute a new 
surface gravity, log g N LTE-o P t, metallicity, [Fe/H] N LTE-o P t, 
and microturbulence, UNLTE-Opt- This was done using the 
same iterative techniques as the LTE-Fe method, that is the 
ionisation balance of the abundance of iron from Fel and 
Fen lines. 

There are, however, three important differences. First, 
the stellar effective temperature was held fixed to the op- 
timal value, Topt. Second, we restricted the analysis to 
Fe lines with excitation potentials above 2 eV, since these 
lines are less sensitive to 3D effects as compared to the 
low-ex citation lines (see the discussion in lBergemann et al.l 
l2012bf ). Third, the abundance of iron from each Fel line 
was adjusted according to the NLTE correction for that 
line at the stellar parameters of the current iteration in the 
procedure. The NLTE corre ctions were determined using 
the N LTE grid computed in iLind. Bergemann. fe Asplundl 
(2012) and applied on-the-fly via a wrapper program 
to MODG. Note that the NLTE cal culations presented in 
ILind. Bergemann. fe Asplundl (|2012T > were analogously cali- 
brated using the ionisation equilibria of a handful of well- 
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known stars. Our extended sample, including more stars 
with direct measurements of surface gravity and effective 
temperature, provide support for the realism of this calibra- 
tion. The grid extends down to log g = 1. We imposed 
a routine which linearly extrapolated the NLTE correc- 
tions to below this value. The results of e xtrapolations were 
checked against NLTE grids presented in lBergemann et al.l 
d2012ah and no significant differences were found. Fur- 
ther. iLind. Bergemann. fc Astalundl (|2012T l found very small 
NLTE corrections for Fen lines. We therefore do not apply 
any correction to the Fen lines. 

Iterations continued until the difference between the av- 
erage abundance of iron from the Fell lines and the NLTE- 
adjusted Fel lines were in agreement (within ±0.05 dex) and 
the slope of the relationship between the reduced EW of the 
Fel lines and their NLTE-adjusted iron abundance was min- 
imised. Sections [5] and [6] describe our final stellar parameter 
estimates for this method. 



4 INITIAL LTE-FE PARAMETERS 

The initial LTE-Fe stellar parameters for our sample stars 
are listed in Table [1] Residuals in the minimizations of this 
technique gave typical internal errors of 0.1 dex in both 
log 3lte-Fc and [Fe/H] L TE-Fc and ~ 55 K in Tlte-Fc As 
we show in the following sections, these small internal er- 
rors can be quite misleading as they are not representative 
of the actual accuracy of stellar parameter estimates. Often, 
especially in metal-poor stars, estimates of T e g, log g, and 
[Fe/H], that result from this method are far too low when 
compared to other, more accurate data (cf., Rll). 



EFFECTIVE TEMPERATURE 
OPTIMISATION 



It was found in 



Bergemann et all 

iLind. Bergemann. fc Asplundl (|2012h 



that 



|2012bl ) 

taking 



and 
into 



account NLTE in the solution of excitation equilibrium 
does not lead to a significant improvement of the stellar 
effective temperature. This was also supported by our test 
calculations for a sub-sample of stars. Fel lines formed in 
LTE or NLTE are still affected by convective surface inho- 
mogeneities and overall different mean temperature/density 
stratifications, whic h are most prominent in strong low- 
excitation Fel lines ( Shchuki na. Truiillo Bueno. fc Asplundl 
120051 : iBergemann et al.ll2012bh . Using ID hydrostatic mod- 
els with either LTE or NLTE radiative transfer thus leads 
to effective temperature estimates that are too low when 
the excitation balance of Fel lines is used (see below). It is 
therefore important that the stellar effective temperature 
be estimated by other means. 



5.1 Three Effective Temperature Scales 

We used three different methods to compute the effective 
temperature. 

The first estimate, Trial, was derived from the wings 
of the Balmer lines, which is among the most reliable 
methods available for th e effective temperature determina- 
tion of FGK stars (e.g., iFuhrmann. Axer. fc Gehrenlll993l ; 




6555 6560 6565 

Wavelength (A) 



6570 



Figure 1. Example fit to H a in the spectrum of the metal-poor 

([Fe/H] 1.0) dwarf, J142911.4-053131. The solid, red curve 

shows the best fit to the data (at 5700 K), while the blue, dashed 
curves represent ±100 K around the fit. 



FuhrmannI Il998l; iBarklem et all |2002|; ICowlev fc Castellil 



20021 ; iGehren et al.ll2006l ; lMashonkina et al.ll200Sft . The only 
restriction of this method is that for stars cooler than 4500 
K, the wings of Hi lines become too weak to allow reliable 
determination of T e ff. Profile fits of H Q and H/j lines were 
performed by careful visual inspection of different portions 
of the observed spectrum in the near and far wings of the 
Balmer lines which were free of contaminant stellar lines. 
Figures [TJ and [5] show two example fits to H Q . Note that 
the Balmer lines were self-contained within a single order 
in each spectrum. Therefore, we did not use neighbouring 
orders for the continuum normalisation. 

Theoretical profiles were computed usin g the SIU 
code with MAFAGS -QDF model atmospheres ( FuhrmannI 
1 19981 ; iGrumil l2004al L Same as ATLAS-0DF (Section [3]), the 
MAFAGS models were computed with Kurucz opacity dis- 
tribution functions, thus the differences between the model 
atmosphere stratifications are expected to be minimal in our 
range of st ellar parameters. Fo r self-broadening of H lines, 
we used the lAli fc Grieml (|l965l ) theory. As shown bv lGruppl 
(|2004al 'l this method successfully reproduces the Balmer line 
spectrum of the Sun within 20 K, and provides accurate stel- 
lar p arameters that agree very well with Hipparcos astrom- 
etry l|Grupp|l2004bf ). The errors are obtained directly from 
profile fitting, and they are largely internal, ±50 to 100 K. 

A key advantage of the Balmer lines is that they are 
insensitive to interstellar reddening, which affects photo- 
metric techniques (see below). However, the Balmer line ef- 
fective temperature scale could be affected by systematic 
biases, caused by the physical limitation of the models. 
The influence of deviations from LTE in the Hi line for- 
m ation in application to co ol metal-poor stars was studied 
bv lMashonkina et al.l j200Sfi. Comparing our results to the 
NLTE estimates by IMashonkina et al.l (|2008T ) for the stars 
in common, we obtain: AT eff (our - M08) = -70 K (HD 
122563 metal-poor giant), AT cff (our - M08) = 50 K (HD 
84937, metal-poor turn-off). The difference is clearly within 
the T fT uncertainties. On the other side, it should be kept 
in mind that the atomic data for NLTE calculations for hy- 
drogen are of insufficient quality and, at present, do not al- 
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Table 1. Atmospheric Stellar Parameter Data 



LTE-Fe 



NLTE-Opt 



Star 






log 9 


[Fe/H] 


v E(B-V) 


TlRPM °T IFtFM 


2rii 


^Bal 




Toff 


a T etf 


log 9 


[Fe/H] 


V 




(K) 


(K) 


(±0.1) 


(±0.1) 




(K) (K) 


(±140 K) 


(K) 


(K) 


(K) 


(K) 


(±0.1) 


(±0.1) 




C0023306-163143 


5128 


58 


2.40 


-2.63 


1.3 




5528 


5400 


100 


5443 


101 


3.20 


-2.29 


0.9 


C0315358-094743 


4628 


40 


1.51 


-1.40 


1.5 




4722 


4800 


100 


4774 


89 


2.06 


-1.31 


1.6 


C0408404-462531 


4466 


40 


0.25 


-2.25 


2.2 




4600 






4600 


90 


1.03 


-2.10 


2.1 


C0549576-334007 


5151 


50 


2.53 


-1.94 


1.3 




5379 


5400 


100 


5393 


82 


3.16 


-1.70 


1.1 


C1141088-453528 


4439 


10 


0.39 


-2.42 


2.1 




4592 


4500 


200 


4562 


123 


1.10 


-2.28 


1.9 



This table is published in its entirety in the electronic edition of the MNRAS. A portion is shown here for guidance regarding its form 

and content. 




Table 2. Effective temperatures determined from direct in- 
tcrfcrometric measurements of angular diameters. 



6555 6560 6565 
Wavelength (A) 



6570 



Figure 2. Example fit to H a in the spectrum of the metal-poor 
([Fe/H] ~ -2.3) giant, J230222. 8-683323. The solid, red curve 
shows the best fit to the data (at 5200 K), while the blue, dashed 
curves represent ±100 K around the fit. 



low accurate quantitative as sessment of NLT E effects in H, 
as elaborately discussed by iBarkleml (|2007l ). Likewise, the 
influe nce of granulation is difficult to assess. iLudwig et al.l 
(2009) presented 3D effective temperature corrections for 
Balmer lines for a few points on the HRD, for which 3D 
radiative-hydrodynamics simulations of stellar convection 
are available. For the SurjJ they find AT cff (3D - ID) « 35 
K, and for a typical metal-poor subdwarf with [Fe/H] = —2, 
AT eff (3D - ID) of the order 50 to 80 K (average over K a , 
H^, and H 7 ). However, in the absence of consistent 3D 
NLTE calculations, it is not possible to tell whether 3D and 
NLTE effects will amplify or cancel for FGK stars. Thus, 
we do not apply any theoretical corrections to our Balmer 
effective temperatures. 

Currently, the only way to understand whether our 
Balmer T e s scale is affected by systematics is by compar- 
ing with independent methods, in particular interferometry. 
We, therefore, computed the Balmer T e g for several nearby 
stars with direct and indirect interferometric angular diam- 
eter measurements. The results are listed in Tabled while 
we plot the difference between our Balmer estimate and that 
from interferometry in Figure [3] Both T e a scales show an 



3 Here, we use their results obtained with «mlt = 0-5 consistent 
with the MAFAGS-0DF model atmospheres adopted here. 



HD 


log 9 


[Fe/H] 


Tint 


°T Int 


^Bal 




Ref 








(K) 


(K) 


(K) 


(K) 




6582 


4.50 


-0.70 


5343 


18 


5295 


100 


a 


10700 


4.50 


-0.50 


5376 


22 


5320 


100 


a 


22049 


4.50 


0.00 


5107 


21 


5050 


100 


a 


22879 


4.23 


-0.86 


5786 


16 


5800 


100 


b* 


27697 


2.70 


0.00 


4897 


65 


4900 


100 


c 


28305 


2.00 


0.00 


4843 


62 


4800 


100 


c 


29139 


1.22 


-0.22 


3871 


48 


4000 


200 


c 


49933 


4.21 


-0.42 


6635 


18 


6530 


100 


b 


61421 


3.90 


-0.10 


6555 


17 


6500 


100 


a 


62509 


2.88 


0.12 


4858 


60 


4870 


100 


c 


84937 


4.00 


-2.00 


6275 


17 


6315 


100 


b* 


85503 


2.50 


0.30 


4433 


51 


4450 


100 


b 


100407 


2.87 


-0.04 


5044 


33 


5025 


100 


b 


102870 


4.00 


0.20 


6062 


20 


6075 


100 


a 


121370 


4.00 


0.20 


5964 


18 


5975 


100 


a 


122563 


1.65 


-2.50 


4598 


42 


4650 


100 


d 


124897 


1.60 


-0.54 


4226 


53 


4240 


200 


c 


140283 


3.70 


-2.50 


5720 


29 


5775 


100 


b* 


140573 


2.00 


0.00 


4558 


56 


4610 


100 


c 


150680 


4.00 


0.00 


5728 


24 


5795 


100 


a 


161797 


4.00 


0.20 


5540 


27 


5550 


100 


a 


215665 


2.25 


0.12 


4699 


71 


4800 


100 


c 


References for Tj nt 


a -ICavrel et al.l 


(2011); 


b - Gaia calibration 



stars from U. Heiter (2012, priv. comm.) (Those marked with a 
[*] have angular diameters determined using the 
su rface-brightness relations f rom iKervella et al j2004r); c - 
iMozurkewich et al.l J2003h ; d - ICreevev et al.l J2012T) 



agreement of 3 ± 60 K for stars with [Fe/H] > — 1, while the 
Balmer estimate is ~ 50 K warmer than Tint at the lowest 
metallicities. These differences are well within the combined 
errors in the interferometric and Balmer measurements. This 
suggests that deviations from ID HE and LTE are either 
minimal, or affect both interferometric and Balmer T e ft in 
exactly same way. Also note, for the stars in common with 
ICavrel et all (|201ll ), our estimates are fully consistent. 

For the second method, we utilised the Tycho-2 and 
2MASS photometry of each star to compute effective tem- 
perature estimates, Tirf m, using the infra r ed flu x method 
(IRFM), as presented in Casagra nde et all l|2010l . hereafter 
C10). Note that in C10 the IRFM calculations were ap- 
plied only to dwarfs and subgiants and validated by com- 
parison with a large body of interferometric angular diame- 
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Figure 3. Comparison of the Balmer effective temperature esti- 
mate to that from interferometry, as listed in Table[2] The Balmer 
effective temperature estimates agree with the interferometric es- 
timate to within ~ 100 K. 



ters. However, the same code can be safely applied to lower 
surface gravities, as shown by compari son with newly de - 
termined angular diameters for giants (|Huber et al.ll2012Ti . 
The advantage of this method is that it is much less sen- 
sitive to model assumptions that are required for spec- 
troscopic analyses. However, the quality of the photomet- 
ric data used to compute IRFM effective temperatures, as 
well as interstellar reddening, can still largely affect the re- 
sult. In our case reddening has been estimated usin g the 
iDrimmel. Cabrera-Layers, fc Lopez-Corredoiral (|2003T ) map, 
with distances derived from our spectroscopic log g. Typical 
values of E(B-V) are around 0.05 mag, although for some 
of the brightest giants the value can be considerably larger 
(see Tabled]). 

Finally, we chose the effective temperature estimates 
from Rll (denoted as Tru) as our third effective temper- 
ature scale. These estimates were based on the Rll cal- 
ibration, which was derived from the the trend between 
[Fe/H]i/rE-Fe and the difference between Tlte-fc and the 
2MASS J-Ks photometric effective temperature for several 
globular cluster stars, Hipparcos stars, and low-reddened 
stars in the Rll sample. In principle, this method should 
yield similar effective temperature estimates to that of 
IRFM, since it util ises the colour-temperature tran s forma - 
tions presented in IGonzalez Hernandez fc Bonifacio! l|2009f ). 
which were based upon their IRFM calculations. The ad- 
vantage is that the Rll calibration relies on the J-Ks 
colour, which is less sensitive to reddening (E(J-Ks) ~ 



0.5E(B-V)). However, note that the J-Ks colour correlates 
only mildly with T e g, and thus calibrations involving that 
index exhibit a rather larger internal dispersion (in our case 
139 K for dwarfs and 94 K for giants) when compared to 
IRFM effective temperature estimates for standard stars. 
Further, effective temperatures computed using the calibra- 
ti on in C10 are typically ~ 40 K hott e r than those computed 
in IGonzalez Hernandez fc Bonifacio! (j200ST ) . A detailed ex- 
planation for this discrepancy is given in C10. 



5.2 Comparisons 

We next applied each of the above methods to our sample 
stars, the values of which can be found in Table[TJ Note that 
for several stars, the Balmer lines fell in the middle of the 
order of the spectrum. The continuum cannot be determined 
with sufficient accuracy in such regions. We therefore did 
not measure the Balmer lines for those stars. Further, not 
all stars in our sample have Tycho-2 photometry estimates. 
We were unable to compute an IRFM effective temperature 
estimate for these stars. 

Figure 0] shows the comparisons between the three ef- 
fective temperature estimates when applied to our sample 
stars. The estimates from Tsal and Tru show remarkable 
agreement, with a difference of only 1 ± 79 K. The IRFM 
effective temperatures, however, are systematically higher 
than both T Ba i and Tru (by 119 ± 215 K and 113 ± 186 K, 
respectively), with an increasing dispersion towards hotter 
effective temperatures. 

It is possible that inherent NLTE and 3D effects could 
be influencing the Balmer effective temperature scale, how- 
ever, we see excellent agreement with Tru and interferomet- 
ric measurements. As stated previously, the effective temper- 
ature estimates compute d in C10 are about 40 K warmer 
than those computed in IGonzalez Hernandez fc Bonifaciol 
(2009). Further, the J-Ks calibration has a large internal 
dispersion. However, the difference between Tirpm and Tru 
extends well beyond these limits. 

It is possible that the uncertainty in the interstellar red- 
dening may be systematically affecting the IRFM estimates. 
In order to test the accuracy of the estimates of redden- 
ing fro m the lDrimmel. Cabrera-Lavers. fc Lopez-Corredoiral 
l|2003h map, we also tried to measure E(B — V) using the 
interstellar Na D lines. However, the majority of the stars in 
our sample had multi-component interstellar NaD features, 
or the feature was not discernible from the stellar Na lines. 
Note, for several of the stars in which we could measure 
single-component interstellar NaD lines, the Drimmel et al. 
estimates and the NaD estimates were on average different 
by < 0.02 mag, which translates to a difference of < 100 K 
in effective temperature. Given these differences, reddening 
alone cannot account for the very large differences (^ 300 K) 
for many stars. 

Instead, the large scatter mostly likely arises from the 
poor quality of the Tycho-2 magnitudes for stars fainter than 
VTycho ~ 9- For the kind of stars analysed in this work, 
B and V magnitudes are the dominant contributors to the 
bolometric flux, as compared to the infrared 2MASS magni- 
tudes. Should a star be matched with a brighter (dimmer) 
source in B and V , then the bolometric flux will be over- 
estimated (underestimated) by a very large amount, and 
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Figure 4. Comparison of the three different T e g scales. Both 
the estimate from Rlf (Tru) and the Balmer lines (7|j a i) agree 
within ~ 80 K, while that from IRFM (Tirfm) is systematically 
higher than both, and shows an increasing dispersion with in- 
creasing effective temperature. Error bars, which show the mean 
error in each T e jj scale, are displayed in the upper left corner. 



IRFM will return a systematically higher (lower) Tirfm es- 
timate. 

Using the sample in C10, it is possible to com- 
pute T e fr using both the Tycho-2 VT yc ho-K a and Johnson- 
Cousins Vjc-Ks calibrations. We plot the difference between 
T(V Ty cho-K s ) and T(Vj C -K s ) as a function of the Tycho-2 V- 
magnitude for the C10 sample (black points) in Figure [5] In 
addition, we have over-plotted the difference between Tirfm 
and Teal for the stars in our present sample (red triangles). 
As shown in the figure, both samples exhibit a large scat- 
ter in the difference in effective temperature estimates at 
VTycho *5 9. Further, the hotter stars are on average among 
the faintest stars in our sample and thus have larger errors 
in VTycho- This is a clear indicator that the poor quality 
of the Tycho-2 photometric measurements for our stars is 
responsible for the discrepancy between Tirfm and Teal- 



5.3 Final Effective Temperature Estimates 

From the comparisons above, Balmer line measurements 
provide the most reliable effective temperature estimates for 
all stars in our sample. In addition, Tru exhibits small dif- 
ferences with respect to Tbsi across all effective tempera- 
tures, as illustrated in Figure [4] In contrast, IRFM effec- 
tive temperatures appear to show a large dispersion, which 




Tycho 



Figure 5. An illustration of the difference between different ef- 
fective temperature estimates as a function of Vxycho • The black, 
filled circles represent the difference between T(Vrycho--Ks) and 
T(Vjq-K s ) for the sample of stars in C10. The red, open triangles 
show the difference between Tirfm an d TR a i for the stars in our 
sample. For Vxycho ^ 9, both samples display an increased scat- 
ter in AT ef f . This clearly indicates that the uncertainty in Vrycho 
is poorly constrained for those stars. 



we attribute to large errors in Tycho-2 photometry, as well 
as uncertainty in the reddening. We therefore adopted the 
mean of only Teal and Tru, weighted according to the in- 
ternal errors from each method, as our final To p t estimate. 
This also serves to reduce the internal error on the final 
optimal effective temperature estimate, which was typically 
< 100 K. 

As noted previously, we could not measure the Balmer 
lines in several of the stars in our sample. For those stars, we 
adopted the Tru estimate. The authors of Rll adopted an 
error of 140 K in their Tru estimate, which was derived from 
the residuals in their calibration. However, the comparison 
between Teai and Tru in Figure [4] suggests that this value 
was overestimated. The mean error in To p t, for those stars 
with both a Teal and Tru estimate, was 90 K. We therefore 
adopted this value for stars with only a single Tru estimate. 
Our final To pt values and corresponding errors can be found 
in Table □ 



6 SURFACE GRAVITY AND [FE/H] IN NLTE 

Using the final values of To p t described in the previous sec- 
tion, we derived the remaining stellar parameters using the 
NLTE-Opt method described in Section [3] We first vali- 
dated this methodology by applying our analysis to a sam- 
ple of 18 "standard" stars, which have Hipparcos parallaxes. 
T he spectra for th ese stars were obtained for the analysis 
m iFulbriehtl l|2000h . and are of similar quality to our sam- 
ple. Both the LTE-Fe and NLTE-Opt atmospheric parame- 
ters for each standard star are given in Table [3] Using the 
Hipparcos parallax and an estimate of the bolometric cor- 
re ction derived from the bolometric flux rel ations presented 
in iGonzalez Hernandez fc Bonifaciol (|2009 ). we also com- 
puted an "astrometric surface gravity" (g n — 4-KGMaT 4 / L) 
for each star, which is listed as log g w in Table [3] Note, 
we computed an astrometric surface gravity using other 
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various flux relations onso. Arribas. fe Martinez-Rogerl 
1 19951 ; ICasagrande et alj|20ld : lTorresll20ld ), finding results 
within ,$0.1 dex of that compu t ed us ing the relation in 
iGonzalez Hernandez &: Bonifacio! (|2009h . Further, we as- 
sumed a mass of 0.8 M© for [Fe/H] < -1 and 0.9 M e 
for [Fe/H] ^ — 1. However, a difference of 0.1 Mq will 
only change the astrometric surface gravity by ~ 0.05 dex. 
The NLTE-Opt surface gravity estimates show a remark- 
able 0.02 ± 0.11 dex agreement with the astrometric sur- 
face gravity, while the LTE-Fe estimates are too low by 
-0.32 ± 0.39 dex. 

Given the agreement, we applied the above analysis to 
our sample stars. The final NLTE-Opt estimates for surface 
gravity and metallicity, as well as for the microturbulence, 
can be found in Table [T] We adopted 0.1 dex error in both 
the surface gravity and metallicity, based on our compar- 
isons with the standard stars above. 



7 NLTE-OPT VS. LTE-FE 

In Figure [6] we compare our final NLTE-Opt stellar pa- 
rameters to those derived using the LTE-Fe method. The 
differences in the estimates of effective temperature, sur- 
face gravity, metallicity, and microturbulence all display 
clear trends with decreasing metallicity. The microturbu- 
lent velocity is underestimated by ~ 0.1-0.2 km s _1 until 
[Fe/H]LTE-Fc ~ —1.8, where ulte-Fo becomes larger than 
f NLTE-Opt - The differences between T e g range from 200 to 
400 K for metal-poor giants, and —50 to —100 K for dwarfs. 
The differences for logy and [Fe/H] reach a factor of 30 
in surface gravity (A logy = 1.5 dex) and a factor of 3 in 
metallicity (A [Fe/H] = 0.5 dex) at [Fe/H] ~ -2.5. 

Figure [7] illustrates how the different LTE-Fe and 
NLTE-Opt results can change the position of each star in 
the log g vs. log(T c ff) plane. In addition, we have included 
sever al evolutionary tracks, computed using the GARSTEC 
code (|Weiss fc S chlattl 2o"o"a). for comparison. Generally, the 
NLTE-Opt estimates of surface gravity and effective temper- 
ature trace the morphology of the theoretical tracks much 
more accurately. Several features are most notable. The 
NLTE-Opt parameters lead to far less stars that lie on or 
above the tip of the red giant branch, and more stars occupy 
the middle or lower portion of the RGB. Also, stars at the 
turn-off and subgiant branch are now more consistent with 
stellar evolution calculations. 

Figures [6] and [7] further prompted us to determine the 
relative importance of the effective temperature scale versus 
the NLTE corrections for gravities and metallicities in the 
NLTE-Opt method. We singled out the effect of the NLTE 
corrections by deriving additional, LTE-Opt surface gravity 
and metallicity estimates using LTE iron abundances com- 
bined with our Topt estimate. Note, as with the NLTE-Opt 
method, Fe lines which have an excitation potential below 
2 eV were excluded. The comparison between these LTE- 
Opt estimates and the final NLTE-Opt estimates are shown 
in Figure [S] As evident from this figure, solving for ionisa- 
tion equilibrium in NLTE also leads to systematic changes 
in the logy and [Fe/H], such that LTE gravities are under- 
estimated by 0.1 to 0.3 dex, whereas the error in metallicity 
is about 0.0 5 to 0.15 dex. These effects ar e cons istent with 
that seen in lLind. Bergemann. fc Asplundl (|2012t ). 
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Figure 6. Comparison of stellar parameters derived using 
the LTE-Fe method and the NLTE-Opt stellar parameters vs. 
[Fc/HJlte-Fc- The difference in effective temperature, surface 
gravity, and metallicity show a large systematic increase with de- 
creasing metallicity. The dual trends seen in A T e g , A log g, and 
A [Fe/H] are a result of the Rll effective temperature calibration, 
in which the authors found that stars with effective temperatures 
less than 4500 K only required a small correction to TLXE-Fe- 
Therefore, these stars stand out in the plots. 
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Figure 7. Stellar positions in the log g vs. T c ff plane in LTE- 
Fe (left panel) and NLTE-Opt (right panel). The curves shown 
are evolutionary track s computed using the GARSTEC code 
llWeiss fc Sc hlattl 200i|). Each track was computed assuming a 
mass of 0.8 Mq and an [Fe/H] of -0.5 (solid, red), -1.0 (long- 
dashed, green), -1.5 (short-dashed, blue), and -2.5 (dot-dashed, 
magenta). Note, those stars around log g = 2.0 that lie away 
from the tracks are most likely horizontal branch stars. 
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Table 3. Hipparcos Star Surface Gravity Comparisons. 
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Figure 8. Comparison of LTE-Opt parameters to NLTE-Opt 
parameters. This shows the underlying effects of the NLTE cor- 
rections. Typically, the LTE-Opt surface gravity and metallicity 
estimates are systematic underestimated by about 0.1 dex and 
0.05 dex, respectively. However, they become increasingly more 
underestimated to about 0.15 dex and 0.3 dex, respectively, at 
the lowest metallicities. 



We thus conclude that reliable effective temperatures 
are necessary to avoid substantial biases in a spectroscopic 
determination of \ogg and [Fe/H], such as displayed in Fig- 
ure [6] We have shown here that, at present, excitation bal- 
ance of Fel lines with ID hydrostatic model atmospheres 
in LTE does not provide the correct effective temperature 
scale, supporting the results by Bergemann et al. (2012). 
On the contrary, Balmer lines provide such a scale. Further- 
more, NLTE effects on ionisation balance are necessary to 
eliminate the discrepancy between Fel and Fell lines, an ef- 
fect that is present, regardless of the adopted T e g. Only in 
this way is it possible to determine accurate surface gravity 
and metallicity from Fe lines. 



8 CONCLUSION 

In this work, we explore several available methods to de- 
termine effective temperature, surface gravity, and metal- 
licity for late-type stars. The methods include excitation 
and ionization balance of Fe lines in LTE and NLTE, semi- 
empirically calibrated photometry (Rll), and the Infra- Red 
flux method (IRFM). Applying these methods to the large 
set of high-resolution spectra of metal-poor FGK stars se- 
lected from the RAVE survey, we then devise a new ef- 
ficient strategy which provides robust estimates of their 
atmospheric parameters. The principal components of our 
method are (i) Balmer lines to determine effective tempera- 
tures, (ii) NLTE ionization balance of Fe to determine log g 
and [Fe/H], and (iii) restriction of the Fel lines to that with 
the lower level excitation potenti al greater than 2 eV to mi n- 
imize the influence of 3D effects (|Bergemann et aT]|2012bl ). 

A comparison of the new NLTE-Opt stellar parame- 
ters to that obtained from the widely-used method of ID 
LTE excitation-ionization of Fe, LTE-Fe, reveals significant 
systematic biases in the latter. The difference between the 
NLTE-Opt and LTE-Fe parameters systematically increase 
with decreasing metallicity, and can be quite large for the 
metal-poor stars: from 200 to 400 K in T c g , 0.5 to 1.5 dex in 
log g, and 0.1 to 0.5 dex in [Fe/H]. These systematic trends 
are largely influenced by the difference in the estimate of 
the stellar effective temperature, and thus, a reliable effec- 
tive temperature scale, such as the Balmer scale, is of critical 
importance in any spectral parameter analysis. However, a 
disparity between the abundance of iron from Fe I and Fe II 
lines still remains. It is therefore necessary to include the 
NLTE effects in Fel lines to eliminate this discrepancy. 

The implications of the very large differences between 
the NLTE-Opt and LTE-Fe estimates of atmospheric pa- 
rameters extend beyond that of just the characterisation 
of stars by their surface parameters and abundance anal- 
yses. Spectroscopically derived parameters are often used to 
derive other fundamental stellar parameters such as mass, 
age and distance through comparison to stellar evolution 
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models. The placement of a star along a given model will 
be largely influenced by the method used to determine 
the stellar parameters. For example, distance scales will 
change, which could affect the abundance gradients mea- 
sured in the Milky Way (e.g., Rfl), as well as the con- 
trove rsial identification of different component s in the MW 
halo |Schonrich. Asplund. fe Casagrandell201ll ; iBeers et al.l 
2012). We explore this in greater detail in the nex t paper of 
this series l|Serenelli. Bergemann. fe Ruchtill2012h . 



ACKNOWLEDGEMENTS 

We acknowledge valuable discussions with Martin Asplund, 
and are indebted to Ulrike Heiter for kindly providing inter- 
ferometric temperatures for several Gaia calibration stars. 
We also acknowledge the staff members of Siding Spring Ob- 
servatory, La Silla Observatory, Apache Point Observatory, 
and Las Campanas Observatory for their assistance in mak- 
ing the observations for this project possible. Greg Ruchti 
acknowledges support through grants from ESF EuroGene- 
sis and Max Planck Society for the FirstStars collaboration. 
Aldo Serenelli is partially supported by the European Union 
International Reintegration Grant PIRG-GA-2009-247732, 
the MICINN grant AYA2011-24704, by the ESF EURO- 
CORES Program EuroGENESIS (MICINN grant EUI2009- 
04170), by SGR grants of the Generalitat de Catalunya and 
by the EU-FEDER funds. 



REFERENCES 

Ali, A. W., & Griem, H. R. 1965, Physical Review, 140, 
1044 

Allende Prieto C, Beers T. C, Wilhelm R., Newberg H. J., 
Rockosi C. M., Yanny B., Lee Y. S., 2006, ApJ, 636, 804 

Alonso A., Arribas S., Martinez-Roger C, 1995, A&A, 297, 
197 

An D., et al., 2009, ApJ, 707, L64 

Arnadottir A. S., Feltzing S., Lundstrom I., 2010, A&A, 
521, A40 

Asplund M., 2005, ARA&A, 43, 481 

Bai G. S., Zhao C, Chen Y. Q., Shi J. R., Klochkova V. G., 
Panchuk V. E., Qiu H. M., Zhang H. W., 2004, A&A, 425, 
671 

Barden S. C, et al., 2008, SPIE, 7014,149 
Barklem P. S., 2007, A&A, 466, 327 

Barklem P. S., Piskunov N., O'Mara B. J., 2000a, A&A, 
355, L5 

Barklem P. S., Piskunov N., O'Mara B. J., 2000b, A&A, 
363, 1091 

Barklem P. S., Stempels H. C, Allende Prieto C, 
Kochukhov O. P., Piskunov N., O'Mara B. J., 2002, A&A, 
385, 951 

Barklem P. S., Stempels H. C, Kochukhov O., Piskunov 
N., O'Mara B. J., 2003, in 12th Cambridge Workshop on 
Cool Stars, Stellar Systems, and the Sun, The Future of 
Cool-Star Astrophysics, ed. A. Brown, G.M. Harper, and 
T.R. Ayres, (University of Colorado), 1103 

Beers T. C, et al., 2012, ApJ, 746, 34 

Bergemann M., Kudritzki R.-P., Plez B., Davies B., Lind 
K., Gazak Z., 2012a, ApJ, 751, 156 



Bergemann M., Lind K., Collet R., Magic Z., Asplund M., 

2012b, arXiv, llarXiv:1207.2455l 
Bruntt H., et al., 2012, MNRAS, 423, 122 
Carollo D., et al., 2010, ApJ, 712, 692 
Casagrande L., Ramirez I., Melendez J., Bessell M., As- 
plund M., 2010, A&A, 512, A54 
Casagrande L., Schonrich R., Asplund M., Cassisi S., 

Ramirez I., Melendez J., Bensby T., Feltzing S., 2011, 

A&A, 530, A138 
Cayrel R., van't Veer-Menneret C, Allard N. F., Stehle C, 

2011, A&A, 531, A83 
Chiavassa A., Freytag B., Masseron T., Plez B., 2011, 

A&A, 535, A22 
Cowley C. R., Castelli F., 2002, A&A, 387, 595 
Creevey O. L., et al., 2012, A&A, 545, A17 
Drimmel R., Cabrera-Lavers A., Lopez-Corredoira M., 

2003, A&A, 409, 205 
Francois P., et al., 2003, A&A, 403, 1105 
Fuhrmann K., 1998, A&A, 338, 161 

Fuhrmann K., Axer M., Gehren T., 1993, A&A, 271, 451 
Fulbright J. P., 2000, AJ, 120, 1841 

Gehren T., Shi J. R., Zhang H. W., Zhao G., Korn A. J., 

2006, A&A, 451, 1065 
Gilmore G., et al., 2012, Msngr, 147, 25 
Gonzalez Hernandez J. I., Bonifacio P., 2009, A&A, 497, 

497 

Grupp F., 2004a, A&A, 420, 289 
Grupp F., 2004b, A&A, 426, 309 

Holmberg J., Nordstrom B., Andersen J., 2007, A&A, 475, 
519 

Huber D., et al., 2012, ApJ. larXiv:1210. 00121 

Ivans 1. 1., Kraft R. P., Sneden C, Smith G. H., Rich R. M., 
Shetrone M., 2001, AJ, 122, 1438 

Johnson J. A., 2002, ApJS, 139, 219 

Kervella P., Thevenin F., Di Folco E., Segransan D., 2004, 

A&A, 426, 297 
Lai D. K., Bolte M., Johnson J. A., Lucatello S., Heger A., 

Woosley S. E., 2008, ApJ, 681, 1524 
Lind K., Bergemann M., Asplund M., 2012, arXiv, 

larXiv: 1207.24541 

Ludwig H.-G., Behara N. T., Steffen M., Bonifacio P., 2009, 

A&A, 502, LI 
Mashonkina L., et al., 2008, A&A, 478, 529 
Mashonkina L., Gehren T., Shi J.-R., Korn A. J., Grupp 

F., 2011, A&A, 528, A87 
McWilliam A., Preston G. W., Sneden C, Searle L., 1995, 

AJ, 109, 275 
Mozurkewich D., et al., 2003, AJ, 126, 2502 
Prochaska J. X., Naumov S. O., Carney B. W., McWilliam 

A., Wolfe A. M., 2000, AJ, 120, 2513 
Ruchti G. R., et al., 2010, ApJ, 721, L92 
Ruchti G. R., et al., 2011, ApJ, 737, 9 
Schonrich R., Asplund M., Casagrande L., 2011, MNRAS, 

415, 3807 

Serenelli A., Bergemann M., Ruchti G. R., 2012, MNRAS, 
in prep. 

ShchukinaN. G., Trujillo Bueno J., Asplund M., 2005, ApJ, 

618, 939 
Sneden C, 1973, ApJ, 184, 839 

Sousa S. G., Santos N. C, Israelian C, Mayor M., Monteiro 

M. J. P. F. G., 2007, A&A, 469, 783 
Steinmetz M., et al., 2006, AJ, 132, 1645 



© 2012 RAS, MNRAS 000, [TV?? 



10 Ruchti et al. 

Thevenin F., Idiart T. P., 1999, ApJ, 521, 753 
Torres G., 2010, AJ, 140, 1158 
Weiss A., Schlattl H., 2008, Ap&SS, 316, 99 
Yanny B., et al., 2009, AJ, 137, 4377 



© 2012 RAS, MNRAS 000, [[]-?? 



